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ABSTRACT 

We report the discovery of a young massive stellar cluster embedded in an extended 
HII region, invisible at optical wavelengths where the extinction is Ay ~ 28 magni- 
tudes, associated with the IRAS source 16177-5018. J, H and nbK imaging photometry 
combined with the K$ 2MASS data show the presence of sources with infrared excess 
emission at 2.2 /im, concentrated in an area of about one square parsec around a massive 
young stellar object identified as the IRAS source. This object has a near-mid infrared 
spectral index betweem 2.2 and 25 (j,m. a(IR) = d \og(\F\)/d logA =4.78, characteristic 
of compact H II regions, with bolometric luminosity, inferred from the integrated near 
to far-infrared flux density of 2.8 x 1O 5 L0, which corresponds to a ZAMS star of about 
42Mq. From the color-magnitude diagram we were able to classify the majority of the 
cluster members as reddened massive stars earlier than spectral type B5. 

Subject headings: stars : formation - stars: pre-main sequence - infrared : stars - ISM: 
HII regions -ISM: dust, extintion 



1. Introduction 

Young massive stars are formed in dense giant molecular clouds, where the gas density may be 
as high as nn = 10 5 — 10 6 cm -3 and the temperature as low as T = 10 — 20 K. For this reason, the 
star forming regions that eventually lie in the core of the molecular cloud can be completely obscured 
at optical wavelengths by many magnitudes of extinction. On the other hand, the dense molecular 
clouds can be traced by molecular emission, like CS and NH3 lines at radio frequencies, which are 
powerful tools for the determination of the physical conditions in the clouds. While OB stars are 
embedded in a molecular cloud, they dissociate and ionize the gas, forming compact HII regions 
seen at radio wavelengths, and heat the circumstellar dust which eventually radiates most of the 



1 Based on observations made at Laboratorio Nacional de Astrofisica/MCT, Brazil 



- 2 - 



stellar luminosity in the far- infrared (FIR). Therefore, the most prominent tracers of young massive 
stars are also compact, or ultracompact H II (UC H II) regions, which have electron densities 
greater than 10 4 cm" 3 and diameters smaller than 0.1 pc. Wood & Churchwell (1989) showed that 
these ultracompact H II regions have characteristic spectral energy distributions in the far-infrared, 
occupying a well delimited region in the log[F v (60fj,m) / F v (12(im)] xlog[F^ (25 fim) / F u (12 fxm)] space. 
Bronfman et al. (1996) made an extensive survey of the CS (2-1) line emission toward IRAS point 
sources with colors characteristics of ultra-compact H II regions and Roman-Lopes et al. (2002) 
detected ammonia emission in the direction of the strongest sources in that survey. With the advent 
of the large bidimensional near-infrared array detectors, the morphological and photometric studies 
of extremely young galactic stellar clusters were greatly benefited. At near-infrared wavelengths (1 
to 2.5 /um) it is possible to probe deep into the dense dust clouds where star formation is taking 
place. At those wavelengths, very young objects present large infrared excess due the presence of 
warm circunstellar dust. 

In this paper we present the discovery of a young stellar cluster of massive stars and its 
associated HII region, in the direction of the IRAS source 16177-5018. We combined 2MASS 
(Two Micron All Sky Survey) data with new NIR observations from Laboratorio Nacional de 
Astrofisica (LNA), Brazil. This work is a part of a survey aimed to the identification of stellar 
populations in the direction of IRAS sources that have colors characteristics of ultracompact HII 
regions (Wood & Churchwell 1989) and strong CS (2-1) line emission (Bronfman et al. 1996). The 
studied region is part of the RCW 106 complex, located in the southern Galactic plane at a distance 
of 3.8 kpc (Caswell & Haynes 1987); it presents compact far infrared sources (Karnik et al. 2001), 
intense NH3 (J,K) = (1,1) line emission (Roman-Lopes et al. 2002) and continuum radio emission 
at 5 GHz as well as hydrogen recombination lines (Caswell & Haynes 1987). 

2. Observations and data reduction 

The imaging observations were performed in June 2001 with the Near Infrared Camera (CamlV) 
of Laboratorio Nacional de Astrofisica, Brazil, equipped with a Hawaii 1024x1024 pixel HgCdTe 
array detector mounted on the 0.6 m Boiler & Chivens telescope. The plate scale was 0.47 arc- 
sec/pixel and the mean values of the PSF full width at half maximum (MFWHM) were 1.2, 1.5 and 
2 arcsec at the J, H and nbK images, respectively. The observations consisted of 8'x 8' frames in 
the direction of the IRAS source 116177-5028 using J, H and nbK filters. The narrow band K filter 
is centered at 2.24956 fim with a FWHM of 0.06333 fj,m, covering only the continuum emission. 
The total integration time was 540 s for the J and H bands and 1620 s for the nbK filter, resulting 
in a sensitivity at 3a of 17.5, 16.8 and 13.4 magnitudes, respectively. The frames were dithered by 
60 arcsec, in order to remove bad pixels, cosmic rays and to eliminate the presence of objects with 
extended emission on the construction of the sky images. Individual frames were dark subtracted 
and flat-fielded before being combined. We obtained a set of "dome flats" for each filter using an 
illuminated white spot in the telescope dome. Sky frames were generated from a median-filtered 
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set of the flattened frames in each band. To produce a combined image from each dithered set of 
images, we aligned the frames using the IRAF 2 subroutines, added the aligned frames and trimmed 
the resulting image to remove bad pixels at the edges of the frames. We used DAOFIND to locate 
stars 4 a above the local background and added to the DAOFIND list all stellar objects missed by 
this routine but found by visual inspection of each image. Because of source confusion within the 
region, photometry was obtained using the point spread function fitting algorithm ALLSTAR in 
the DAOPHOT package (Stetson 1987). For the J and H images the adopted PSF fitting radius 
were 2.0 and 3.2 pixels respectively and the PSF radius was 11 pixels. For the nbK image the PSF 
fitting radius was 4.2 pixels and the PSF radius was 13 pixels. The local sky was evaluated in an 
annulus with an inner radius of 8 pixels and a width of 20 pixels. The values of the completeness 
limits for J, H and nbK bands are 16, 15.5 and 12.2 magnitudes respectively and were derived from 
the point where the number of detected sources of magnitude m, N(m) deviates from a straight 
line in the log(iV) versus m diagram as can be seen in figure 1. Photometric calibration was done 
by observing standards from the list of Elias (1982), at several air masses before and after each set 
of integrations. Direct equivalence was assumed between our nbK magnitudes and the standard K 
values. 

Photometry from 2MASS All Sky Point Source Catalogue 3 in the J, H and K$ filters became 
available recently. The Kg filter, centered at 2.17 /im, has a bandpass of 0.32 (j,m. The com- 
pleteness limits for this survey are 15.8, 15.1 and 14.3 magnitudes at J, H and K$, respectively 
(Egan et al. 2001). Since the 2MASS K$ band photometry has a completeness limit greater than 
our nbK photometry, we decided to use the 2MASS catalogue astrometry and magnitudes in this 
filter to study the stellar population. Even though the 2MASS survey was conducted with a 1.3 m 
telescope, the sensitivity and resolution of our IR camera compensates our smaller collecting area, 
resulting in more detection in the J and H bands, as it will be shown in section 3.2. For this reason 
we only used the more sensitive 2MASS Ks photometry. 

We also checked the consistence of our photometry comparing it with the data of the 2MASS 
survey. We took all point sources detected by 2MASS in an area of about 50 square arcmin around 
IRAS 16177-5018 and constructed a list of sources common to both surveys in the three filters. 
This procedure gave us a total of 612, 781 and 232 sources in the J, H and K bands respectively. 
We constructed graphs for M.-2 mass versus M.c a miv J ■, H and K magnitudes respectively, which are 
shown in figure 2. We see a good linear relation between the two systems, with a slope of 1 and a 
dispersion that increases with the magnitude. 



2 IRAF is distributed by the National Optical Astronomy Observatories, which is operated by the Assotiation of 
Universities for Research in Astronomy, Inc. under contract to the National Science Foundation 

3 http://www. ipac.caltech.edu/cgi-bin/gator/nph-dd 
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3. Results &: Discussion 

The infrared images, especially in the H , nbK and K$ bands, show the presence of a small 
infrared nebula (~ 70" x 70" square arcsec) at the IRAS coordinate, with an apparent concentration 
of embedded stars, suggesting the presence of a cluster. In figure 3 we show a combined pseudoreal 
colour (J represented in blue, H in green and nbK in red) image of the whole field, and amplified 
images of the nebular region in the J, H and nbiT bands. 

In figure 4 we present the (J — H) versus (H — Kg) diagram for all stars detected in the 
LNA J and H images and in the 2MASS K$ image or only at H and K$, using the completeness 
limit for the J magnitude, together with the position of the main sequence, giant branch and 
reddening vectors for early and late type stars. The majority of the sources have colors of reddened 
photospheres but some of them appear to the right of the reddening line for early type stars, 
showing "excess" at 2.2 /j,m. When we combine the results shown by the color-color diagram with 
the spatial distribution of the sources, we verify that all stars situated to the right of the early type 
reddening vector with H — K$ > 1.0 are located in the nebular region. 

It is well established that very young pre-main sequence objects present large infrared excess 
due to the presence of warm circumstellar dust. Our result suggests that the stellar population 
in the direction of the IRAS source is very young, and from their position in the (J — H) versus 
(H — Ks) diagram they are probably young massive stellar objects. 

In order to estimate the amount of interstellar extinction in the direction of the IRAS source 
we computed the average value of the color excess E{H — K) using the 11 sources from Table 1 
(IRS 1, 2, 7, 8, 11, 18, 22, 23, 26, 28 and 31) that lie along the reddening line for early type stars 
in the color-color diagram. We assumed that the intrinsic (H — K)q colors are those of luminous 
early type stars (-0.05 for 06-08V types as given by Koornneef, 1983). The mean color excess 
obtained was < E(H — K) > = 1.75, which for a standard reddening law (Rieke & Lebofsky 1985) 
corresponds to a median visual extinction < Ay 28. 

3.1. The IRAS source 

In Figure 5 we present a contour map constructed from the LNA nbK image, which shows the 
region around IRAS16177-5018. The IRAS coordinate has an intrinsic error delimited by the ellipse 
plotted in the figure. Any object inside the ellipse could be identified with the IRAS source and 
in fact, more than one of these objects may be contributing to the detected infrared flux . A more 
accurate position for the IR source was obtained from the Midcourse Space Experiment - MSX 
source catalog 4 . The MSX surveyed the entire Galactic plane within | b |< 5° in four mid- infrared 
spectral bands centered at 8.28, 12.13, 14.65 and 21.34 /j,m, with image resolution of 20 arcsec and 



4 http://www. ipac.caltech.edu/ipac/msx/msx. html 
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a global absolute astrometric accuracy of about 1.9 arcsec (Price et al. 2001). We found only one 
MSX source, with coordinates a(J2000) = 16 h 21 m 31.4 s , <5(J2000) = -50 d 25 m 04 s within the IRAS 
uncertainty ellipse; it coincides with the star we labeled IRS7 in Table 1. In fact this is the only 
MSX source in an area of about 9 square arcmin around the IRAS point source with detected flux 
at the four bands. 

In Figure 6 we plotted the near to far-infrared spectral energy distribution of IRS7, without 
any correction for absorption. We can see that the MSX and IRAS flux densities agree very well 
at the common wavelengths (12 and 21-25 /im), indicating that only one star is responsible for the 
IRAS emission. The spectral index in the near and mid-infrared (betweem 2.2-25 /im), defined as 
a = d \og(\F\)/d logA is a = 4.78, characteristic of extremely young stellar objects in a earlier 
stage of evolution, when the star is deeply surrounded by a thick dusty envelope. 

We integrated the observed flux density between 1.25 and 100 /im and, assuming a distance 
of 3.8 kpc, we obtained a luminosity L = 2.8 x 10 5 Lq. Since most of the energy is emitted in 
the infrared, this can be considered a lower limit to the bolometric luminosity of the embedded 
star, reprocessed by a dusty envelope. PMS models (Iben 1965) show that massive stars (above 9 
M ) evolve at constant luminosity from the Hayashi track to the Main-Sequence. Therefore, the 
measured luminosity gives directly a lower limit to the stellar mass of 42 M Q , corresponding to a 
ZAMS star earlier than 05.5 (Hanson et al. 1997), independent of its PMS evolutionary stage. 

3.2. Cluster Population 

In order to examine the assumption that a cluster of embedded young stars is present in 
the core of the molecular cloud, we analyzed the stellar population in two delimited regions; one 
that we labeled "nebulae" that contains the cluster and another that we labeled "control", as 
ilustrated in Figure 7. In that figure we represented the position of all objects detected in the H 
band by crosses. We can see that the small region labeled "cluster region" show a concentration 
of sources. As mentioned in section 3, all the sources with excess emission at 2.2 /im belong to 
the nebular region, which has an area about 1 square arcmin (~ 1 x 1 sq parsec at a distance 
of 3.8 kpc), compatible with the sizes of other stellar clusters studied with near infrared arrays 
(Persi et al. 1994; Tapia et al. 1996; Persi et al. 1997). 

In order to separate the cluster sources not detected at J band from the field stars, we con- 
structed a comparative H — K s color distribution diagram for the control and "cluster" regions 
(with the counts for the control region normalized to the "cluster" region area), shown in figure 8 . 
From this diagram we can also see the existence of an excess of sources with H — K$ > 1.0 in the 
direction of the nebulae. 

Using these results, the following criteria was used to select the "cluster member candidates" 
for the nebulae region: a) sources with H — K$ > 1.0 that lie to the right or on the reddening vector 
for early type stars in the color-color diagram; b) sources with H — K$ > 1.0 not detected at the J 



-6- 



band and c) sources detected only at Ks band. For the b) and c) criteria, the contamination from 
background stars, mainly late-type giants, should be very small. In fact for an average extinction Ay 
of 28 magnitudes, the extinction in the K band Ak is about 3.1 magnitudes (using the interstellar 
reddening law taken from Rieke & Lebofsky, 1985), and for a distance to the cloud of 3.8 kpc, these 
sources would have K > 16. Table 1 shows the coordinates and photometry of all selected sources 
(J and H from both LNA and 2MASS and K s and coordinates from 2MASS). From the 41 stars 
with measured Ks magnitudes we were able to measure the LNA magnitudes in 35 of them at H 
and 22 at J bands, while the 2MASS survey only detected 20 objects at H and 6 at J bands. 

Further information about the nature of the selected objects in Table 1 can be extracted from 
Ks versus (H — Ks) color-magnitude diagram shown in Figure 9. The locus of the main-sequence 
for class V stars it is also plotted, with the position of each spectral type earlier than AO V 
indicated by filled circles. The intrinsic colors were taken from Koornneef (1983) while the absolute 
K magnitudes were calculated from the absolute visual luminosity for ZAMS taken from Hanson et 
al. (1997), for a distance of 3.8 kpc (Caswell & Haynes). The reddening vector for a ZAMS BO V 
star, taken from Rieke & Lebofsky (1985), is shown by the dashed line with the positions of visual 
extinctions Ay = 10, 20 and 30 magnitudes indicated by filled circles. From this diagram we can 
see that the majority of sources can be classified as reddened massive stars earlier than spectral 
type B5, forming an O-B cluster embedded in the molecular cloud. However, the spectral types 
inferred from the color-magnitude diagram are only upper limits when the stars present emission 
"excess" in the NIR. 

We made a crude estimate of the star formation efficiency using the mass of gas measured 
by Karnik et al. (2001). They mapped the RCW 106 region at 150 and 210 /xm and detected a 
compact region in the direction of the IRAS source 16177-5018. They modeled the observed mid 
and far-infrared SED assuming a constant gas to dust relation M gSuS /M^ ust = 100 and from the 
integrated luminosity about 1O 5 L0 they obtained a mass of dust Md ust ~ 30 — 35M resulting in 
a mass of gas M gas « 3000 — 3500 M§. From this value we estimated a lower limit to the star 
formation efficiency (SFE), using the lower limit to the number of cluster members detected as 41 
and an average stellar mass of 15M , which corresponds to a B0 V star (Hanson et al. 1997): 

SFE= ^ Mst Z = 0.17 (1) 

M gas + M stars 

This is a lower limit because the mass of gas assumed for the nebular region may be overestimated 
since it can include a larger volume than our star forming region and also because stars fainter 
than our detection limit are not included. 

We also estimated the number of Lyman-continuum photons produced in the star forming 
region, from the spectral types of the stars that do not show "excess" in the color-color diagram, 
using the relation between spectral type and number of Lyman-continuum photons given by Hanson 
et al. (1997). We obtained 1.6 x 10 50 photons s -1 . 
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3.3. The infrared nebula 

In figure 10 we shown contour maps of the infrared nebulae region at J, H and x\bK bands 
obtained from our infrared images. The contours were calibrated in flux with the values starting at 
0.18, 0.48 and 2.8xl0 -4 Jy/beam respectively. From the nbif contour map we estimated a lower 
limit to the flux density of the emission nebula in the K band by measuring the area between 
contours and multiplying by the value of the corresponding flux density per unit area and we found 
S(K) = 0.62 Jy. Assuming that the emission process is thermal bremsstrahlung from an optically 
thin region, we can compare the IR with the radio flux density at 5 GHz given by Caswell & 
Hynes (1987) to obtain an upper limit to the K-band absorption. Assuming constant density and 
temperature across the cloud and local thermodynamic equilibrium, the flux density S(u) due to 
free-free emission can be written as: 

S(u) = r u B u (T)n (2) 
where t v is the optical depth at frequency u, B V (T) is the Planck function 

B ^ T ) = C 2 exp(WfcT) - 1 (3) 
and £1 is the solid angle of the source given by: 

n = 7t(L/2D) 2 (4) 

where L is the diameter of the ionized cloud and D the distance to the observer. The optical depth 
at a given frequency v is: 

r v = a s L (5) 
where ag is the free- free absorption coefficient (cgs) taken from Rybick (1979): 

3.7 x 10 s [1 - exp(-hv/kT)} (n e n,) g s (u, T) 
as = ^ z -2 T i/2 ( 6 ) 

gs(v,T) is the Gaunt factor obtained from Clayton (1968) 

For two frequencies v\ and the ratio of the corresponding flux densities S{y{) and 5(^2) 
may be calculated from: 

— — = exp[h{v 2 ~ vi) kT\— — (7) 

S{V2) gs{V2,T) 

For a flux density in the 5 GHz continuum of 32.4 Jy and an electron temperature T e = 6300 K, 
as given by Caswell & Haynes (1987) we find an expected flux density at the K band of 11.5 Jy. The 
measured value of 0.62 Jy gives an upper limit to the K band absorption of 3.1 magnitudes. Using 
the interestellar extinction law given by Rieke & Lebofsky (1985) we obtained a visual extinction 
Ay k, 27.7 in very good agreement with the value obtained in section 3.0 from the median (H — K) 
color excess. 
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We observe that is possible that the detected infrared flux at 2.2/im has a contribution from 
light scattered by dust, but taking into account that the ionized region detected at radio wavelengths 
coincides with the near infrared nebulae and the high number of massive stars inferred from the 
color-magnitude diagram, we believe that the major contribution to the K band integrated flux is 
probably due to free- free emission from the ionized region. 

We also obtained the number of ionizing Lyman continuum photons Ni y from the radio con- 
tinuum flux density using the expression derived by Rubin (1968): 



where v is in MHz and /j is the fraction of helium recombination photons that are energetic 
enough to ionize hydrogen (/« « 0.65). We assumed, as in Simpson & Rubin (1990) the ratio 
He/H = 0.1 and the average volume of He + as half the volume of H + . Using Equation 7 we found 
Nhy ~ 5.2 x 10 49 photons s _1 . This value is compatible with the number of Lyman continuum 
photons available from the main sequence stars derived in the previous section. 



Near-IR imaging photometry in the direction of the IRAS source 16177-5018 revealed a young 
massive stellar cluster and its associated HII region embedded in a dense molecular cloud . We 
detected 41 cluster member candidates concentrated in an area of 1 square parsec, but probably 
other less massive sources were missed because the completeness limit. 

The NIR counterpart of the IRAS point source was identified using more accurate positions 
from the MSX catalogue. This source has a near-mid spectral index d \og{\F\)/d logA = 4.78 
indicating that the IRAS source is very young, with a lower limit to the bolometric luminosity 
inferred from the integrated near to far-infrared flux density of about 2.8 x 10 5 L & , which corresponds 
to a 42 M Q star. 

The average absorption found for the embedded early type stars {Ay ~ 28) coincides with the 
absorption of the HII region, derived by comparing the integrated extended flux density at the K 
band to what is expected at this band from the observed radio emission. 

The number of Lyman continuum photons expected from the 11 stars that lie along the red- 
dening line for early-type stars about 1.3 x 10 50 s" 1 , enough to produce the 5 GHz flux density 
and near IR emission, assumed to be originated by free-free emission. 

This work was partially supported by the Brazillian agencies FAPESP and CNPq. We acknowl- 
edge the staff of Laboratorio Nacional de Astrofisica for their efficient support and to F. Jablonsky 
who was the PI of the CamlV project. We also thank the anonymous referee for constructive 
suggestions that improved the quality of this paper. This publication makes use of data products 




5.99 x 10 47 ' S{v)D 2 T-° Ab v 0A 
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(8) 



3.4. 



Conclusions 
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from the Two Micron All Sky Survey, which is a joint project of the University of Massachusets 
and the Infrared Processing and Analysis Center/California Institute of Technology, funded by the 
National Aeronautics and Space Administration and the National Science Foundation. 
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Table 1. List of the selected near-infrared sources 
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Fig. 2. — The magnitudes comparative diagram Mjhk (2MASS) x MjHK(CamIV). The continu- 
ous line shows the relation expected if the two photometric systems were equal. 
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Fig. 3. — Color-color diagram of all sources detected at the three passbands. The locus of the main 
sequence and giants branch are shown by the continuous lines taken from Koornneef (1983), while 
the two parallel lines (dashed and dotted) follow the reddening vectors taken from Rieke & Lebofsky 
(1985). The location (crosses) of Ay = 0, 10, 20, 30 and 40 magnitudes of visual extinction are 
indicated by black bold numbers. We also plotted the sources only detected at H and 2MASS Kg 
images using the completenees limit of the J photometry (about 16 mag). These sources have lower 
limit to the J-H color indicate by arrows. The cluster members candidates taken from table 1 are 
labeled by the red numbers. 
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Fig. 4. — nbK band contour map of the infrared nebulae associated with IRAS16177-5018. The 
contours start at 1.8xl0 -4 Jy/beam, with the same intervals (the beam size is 2 x 2 pixels). The 
positions of selected infrared sources taken from table 1 are indicated by black numbers. It is also 
indicated the location of the IRAS coordinate elipse error (black dotted line). 
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Fig. 5. — Spectral energy distribution of the IRS7 infrared source. The flux for the three near 
infrared bands (filled circles) were taken from our survey. The mid infrared data (filled squares) 
were taken from MSX catalogue (bands A=8.28/mi, C=12.13/xm, D=14.65 / um and E=21.34/xm) 
while the far infrared data (filled down triangles) were taken from IRAS point source catalogue. 
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Fig. 6. — Diagram of the spacial distribution of sources detected at LNA's H band . The two 
comparative regions a (nebulae) and b (control) are delimitated by the vertical continuum line. 
The "cluster region" area is delimited by the dotted square at the up right side. 



-17- 



60 



50 



"Nebulae" 
control 



40 - 

ST 

X 30 - 
z 

20 - 




H-K 



Fig. 7. — Comparative diagram of the distribution of the H — K color for the " control" (continuous 
line) and "nebulae" (dotted line) areas . 
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Fig. 8. — Ks x (H — Ks) color-magnitude diagram of the selected sources in table 1. The locus of 
the main sequence (dwarfs) at 3.8 kpc is shown by the continuous line. The intrinsic colors were 
taken from Koornneef (1983) while the absolute K magnitudes were calculated from the absolute 
visual luminosity for ZAMS taken from Hanson et al. (1997). The reddening vector for a B0 ZAMS 
star (dotted line) was taken from Rieke & Lebofsky (1985); Also are indicated the location (filled 
circles) of Ay = 10, 20 and 30 magnitudes of visual extinction. The sources IRS1-IRS35 are labelled 
by red numbers. 
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Fig. 9. — J, H and nbK contour maps of the infrare nebulae region. The contours start at 1.8xlCT 5 
Jy/beam in J, 4.8xlCT 5 Jy/beam in H and 2.8xlCT 4 Jy/beam in nbK, and are spaced by the 
same fluxes per beam (a beam is 2 x 2 pixels) . 
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